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Abstract 



OO ! The Navier-Stokes systems for compressible fluids with density-dependent 

' viscosities are considered in the present paper. These equations, in particular, 

include the ones which are rigorously derived recently as the Saint- Venant system 
^ , for the motion of shallow water, from the Navier-Stokes system for incompressible 

H I flows with a moving free surface [14]. These compressible systems are degenerate 

when vacuum state appears. We study initial-boundary- value problems for such 
systems for both bounded spatial domains or periodic domains. The dynamics 
of weak solutions and vacuum states are investigated rigorously. 

First, it is proved that the entropy weak solutions for general large initial 
data satisfying finite initial entropy exist globally in time. Next, for more regular 
initial data, there is a global entropy weak solution which is unique and regular 
with well-defined velocity field for short time, and the interface of initial vacuum 
propagates along particle path during this time period. Then, it is shown that 
for any global entropy weak solution, any (possibly existing) vacuum state must 
vanish within finite time. The velocity (even if regular enough and well-defined) 
blows up in finite time as the vacuum states vanish. Furthermore, after the 
vanishing of vacuum states, the global entropy weak solution becomes a strong 
solution and tends to the non-vacuum equilibrium state exponentially in time. 
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Vanishing of vacuum states and blow-up phenomena 



1 Introduction 

The compressible isentropic Navier-Stokes equations, which are the basic models 
describing the evolution of a viscous compressible fluid, read as follows 



where x eQ C R^,t e (0,T), D{u) = (Vu + (Vu)^)/2, and p(p) = ap^,a > 0,7 > 1, 
the viscosity coefficients ^ are assumed to satisfy /i > and ^ + 2/x/iV > 0. 

If n and ^ are both constants, there is huge literature on the studies of the global 
existence and behavior of solutions to (11.11) . For instance, the one-dimensional (ID) 
problems were addressed by Kazhikhov et al [24] for sufficiently smooth data, and by 
Serre [46,47] and Hoff [18] for discontinuous initial data, where the data are uniformly 
away from the vacuum; the multidimensional problems (11.11) were investigated by Mat- 
sumura et al [35-37], who proved global existence of smooth solutions for data close 
to a non- vacuum equilibrium, and later by Hoff for discontinuous initial data [19], and 
more recently, by Danchin [11], who obtained existence and uniqueness of global solu- 
tions in a functional space invariant by the natural scaling of the associated equations; 
and for the existence of solutions for arbitrary data(which may include vacuum states). 
Lions [29-31] (see also Feireisl et al [13]) obtained global existence of weak solutions - 
defined as solutions with finite energy - when the exponent 7 is suitably large, where 
the only restriction on initial data is that the initial energy is finite, so that the density 
is allowed to vanish. 

Despite the important progress, the regularity, uniqueness and behavior of these 
weak solutions remain largely open. As emphasized in many papers related to com- 
pressible fluid dynamics [9,10,18,20,21,24,27,45,46,48,50,52], the possible appearance 
of vacuum is one of the major difficulties when trying to prove global existence and 
strong regularity results. Hoff and SmoUer [21] proved that weak solutions of the com- 
pressible Navier-Stokes equations (11.11) in one space dimension do not exhibit vacuum 
states in a finite time provided that no vacuum is present initially under fairly gen- 
eral conditions on the data. Such a result was extended to the spherically symmetric 
case in [54] recently. On the other hand, the results of Xin [52] showed that there is 
no global smooth solution to Cauchy problem for (11.11) with a nontrivial compactly 
supported initial density, which gives results for finite time blow-up in the presence 
of vacuum. It is also proved in [27] that for bounded domain even one point initial 
vacuum shall cause the global strong solutions to ID (II. ip to blow up as time goes to 
infinity provided the initial data satisfies some compatibility conditions. 

The independence of viscosities on density makes it possible to trace either the 
particle path or the trajectory of vacuum state. This, however, leads to the failure of 
continuous dependence of weak solutions containing vacuum state on initial data [20]. 
For the case that the density changes continuously across the interfaces separating the 
gas and vacuum, the global existence and uniqueness of weak solutions was obtained 



pt + div{pu) = 0, 

{pu)t + div{pu ® u) — 2div{nD{u)) — V{^divu) + Vp{p) = 0, 
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in [33] , where the authors obtained that velocity is smooth enough up to the interfaces 
which are particle paths separating the gas from the vacuum and that the support of 
gas density expands outside as well as interface connecting gas and vacuum moves at 
an algebraic rate. 

Thus, viscous compressible fluids near vacuum should be better modeled by the 
compressible Navier-Stokes equations with density-dependent viscosities, as was de- 
rived in the fluid-dynamical approximation of Boltzmann equation for dilute gases. 
Further, as was first pointed out and investigated by Liu-Xin-Yang in [32], in the 
derivation of the compressible Navier-Stokes equations from the Boltzmann equation 
by the Chapman-Enskog expansions, the viscosity depends on the temperature, which is 
translated into the dependence of the viscosity on the density for isentropic flows. More- 
over, it should be emphasized that a one- dimensional compressible flow model, called 
the viscous Saint- Venant system for laminar shallow water, derived rigorously from 
incompressible Navier-Stokes system with a moving free surface by Gerbeau-Perthame 
recently in [14], has the form: 



with the viscosity coefficients given by yu(p) = ^(p) = ap/3 for a given positive constant 
a. Indeed, such models appear naturally and often in geophysical flows [2,6,7]. 

In the case of one-dimensional problem with p = ^ = ap" for some positive constants 
a and a, the well-posedness of the Cauchy problem has been studied by many authors 
for not only initially compact-supported density but initially non-vacuum density. In- 
deed, for the initially compact-supported density case, the local (in time) well-posedness 
of weak solutions to this problem was first established by Liu-Xin-Yang in [32], where 
the initial density was assumed to be connected to vacuum with discontinuities. This 
property, as shown in [32], can be maintained for some finite time. And the global ex- 
istence of weak solutions, together with a G (0, 1) and the density function connecting 
to vacuum with discontinuities, was considered by many authors, see [23,43,55] and 
the references therein. It is noticed that the above analysis is based on the uniform 
positive lower bound of the density with respect to the construction of the approximate 
solutions. On the other hand, if the density function connects to vacuum continuously, 
there is no positive lower bound for the density and the viscosity coefficient vanishes 
at vacuum. This degeneracy in the viscosity coefficient gives rise to new difficulties in 
analysis because of the less regularizing effect on the solutions. Yang et al [56] first 
obtained a local existence result for this case under the free boundary condition with 
a > 1/2. The authors in [12,51,57] obtained the global existence of weak solutions 
with a e (0,1/2). 

However, almost all above results concern mainly with free boundary problems, 
and for the global existence results, the choices of viscosity do not fit the important 
physical model, the shallow water equation (11.21) with p(p) = p (namely a = 1). For the 
constant viscosity case, one has known not only that the vacuum state will not develop 



Pt + {pu)x = 0, 

{pu)t + {pu^)x - a{pux)x + (p^)x = 0, 
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later on time if there is no vacuum state initially [21], but also that the separate two 
initial vacuum states shall not meet together in a finite time [54], and that one point 
initial vacuum causes strong solutions to blow up [27] at infinity as well. However, little 
is known on the dynamics of the vacuum states of weak solutions to the compressible 
Navier-Stokes equations (11.11) with density-dependent viscosity on bounded domain. 
And in particular, it is not clear yet how the vacuum states evolve with respect to 
time and whether the initial vacuum states shall exist all the time or not for weak 
solutions. The study of these important dynamical problems about vacuum states 
is rather difficult because the nonlinear diffusion is degenerate as vacuum appears, 
which is quite different from the case of constant viscosity. This causes the loss of 
information about the velocity and makes it difficult to trace the evolution of vacuum 
states in general. Across the interface (or vacuum boundary), it is usually difficult to 
obtain enough information about velocity even if considering specific cases such as point 
vacuum or continuous vacuum of one piece. It is important to get enough information 
about the velocity since the fiow particles transport usually along particle path, and 
all interfaces of vacuum, such as free boundaries [23, 33] which can be observed and 
dealt with, also move along the trajectories determined by velocity field. 

For the multidimensional case, Vaigant et al [49] first proved that for the 2D case 
and for the case /i is a constant and ^(p) = ap^, with a > 0, /3 > 3, (II. ip with periodic 
boundary condition has a unique strong and classical solution with density away from 
vacuum. More recently, Bresch and Desjardins [2,3,6,7] (see also [38]) have made 
important progress. Under the condition that ^{p) = 2(p'(p)p — /i(p)), they estabhsh 
a new Bresch- Desjardins (BD) entropy inequality which can not only be applied to 
the vacuum case but also used to get the compactness results for (II. ip which extended 
the compactness results due to Lions [29-31] to the case 7 > 1. On the other hand, 
the constructions of the approximation solutions does not seem routine in the general 
case of appearance of vacuum. However, it should be noted that recently Bresch et 
al [4-6] have made significant progress on the construction of approximate solutions 
and existence of global weak solutions to the mult i- dimensional compressible Navier- 
Stokes equations and the 2D shallow water model in the case that there is either a 
drag friction (for barotropic compressible fiows) or a cold pressure (for viscous and 
heat conducting fiows). In the case that there is neither drag friction nor cold pressure 
included, Guo-Jiu-Xin [16] recently have shown how to construct approximate smooth 
solutions and obtain the global existence of weak solutions via the BD entropy to 
the (2D and 3D) barotropic compressible Navier-Stokes for the spherical symmetric 
initial data. There are also other recent interesting applications of the BD entropy 
to one-dimensional compressible Navier-Stokes equations with degenerate viscosities, 
for instance, on the global existence and long time behavior of weak solutions for free 
boundary problem [15,17], or the existence and uniqueness of global strong solution 
away from vacuum in real line [39]. 

We study mainly the initial-boundary- value problem (IB VP) for (11.11) . where /i = p° 
with a > 1/2, on spatial one-dimensional bounded spatial domains or periodic domains. 
This contains the physical important model for shallow water equations (11.20 . The 
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choice of a > 1/2 is necessary in order to consider the dynamics of vacuum states since 
it allows the existence of initial vacuum in Eulerian coordinates as one can see later. 

We first establish the global existence of entropy weak solutions for the compressible 
Navier-Stokes equations (11 .ip . with pressure p = and 7 > a/2, for general initial 
data with finite entropy and vacuum. The key in our analysis is the construction non- 
vacuum approximate solutions so that we can make use of the stability analysis in [38], 
where the Bresch-Desjardins (BD) entropy inequality (see [2,3,6,7]) was used to obtain 
the compactness results. Our construction of the approximate solutions is strongly 
motivated by the previous work of Jiang-Xin-Zhang [23] about the existence of global 
weak solutions to one-dimensional compressible Navier-Stokes with free boundary as 
vacuum interface. In general, it seems rather difficult to investigate the dynamics of 
vacuum states due to the degeneracy of nonlinear diffusion and the density function 
connecting to vacuum continuously. Therefore, we further consider the cases of more 
regular initial data containing point vacuum or continuous vacuum of one piece, and 
we show that there is a global entropy weak solution which is unique and regular with 
well-defined velocity field at least for short time, and the vacuum states remain for 
the short time. Then, we use some ideas due to [22, 27, 28] to prove that any possible 
vacuum state in such global weak solutions which satisfy the BD entropy must vanish 
within finite time. This shows that such short time structure and vacuum states of 
weak solutions can not be maintained all the time. And as the vacuum states vanish, 
the spatial derivative of velocity (if it exists) has to blow up even if the velocity is 
regular enough and well-defined before. After the vanishing of vacuum states, we can 
redefine the velocity field and recover the nonlinear diffusion term in terms of density 
and velocity. In addition, the global entropy weak solution is shown to become a strong 
solution and tends to the non-vacuum equilibrium state exponentially in time. This 
phenomena, applied to the compressible shallow water equations (II. 2p . seems to be 
never observed for the compressible Navier-Stokes equations before. 



The rest of the paper is as follows. In section [21 the main results about the vanishing 
of vacuum states and blow-up phenomena of global entropy weak solutions for the 
compressible Navier-Stokes equations are stated. The global existence of entropy weak 
solutions for general large initial data with vacuum states allowed is proven in section [3l 
The short time structure of global entropy weak solution with initial one point vacuum 
state or initial continuous vacuum states of one piece are investigated in section HI 
In section |5l we show the vanishing of vacuum states and blow-up phenomena of any 
global entropy weak solution within finite time and analyze the regularity and large time 
asymptotic behavior of global entropy weak solutions after the vanishing of vacuum 
states. 
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Vanishing of vacuum states and blow-up phenomena 



2 Main results 

We consider the initial-boundary-value-problem (IBVP) for the ID compressible 
Navier-Stokes equations with density-dependent viscosity 

Pt + (pu), = 0, (2.1) 
{pu)t + {fM^ + p(p))x - (/i(p)Mx)x = 0, (2.2) 

with p > the density, pu the momentum. The pressure and viscosity are assumed to 
have the form: 

Pi.p) = «lP^7 P{p) = «2P" 

where 7 > 1, ai > 0, a2 > 0, and a > 1/2 are constants, and for simplicity we set 
fli = = 1. 

The initial data is given for the density p and the momentum pu 

p(x, 0) = Pq{x) > 0, pu{x,0) = mo{x), x E fl, (2.3) 

where the domain Q is chosen as unit interval denoting the spatial domain (0, 1) or 
periodic domain with period length one, and throughout the present paper the initial 
data is assumed to satisfy 



po > a.e. in Q, po e L\n), (po"'/'). G L^Q), 
mo = 0, a.e. on{xen\ po{x) =0}, ^ G L\n). 



(2.4) 



Remark 2.1 Note here that the condition (\2.4^ implies 



PoeL^iQ), polog^poeL\Q). (2.5) 

It should be clear that a large class of initial data satisfy the conditions in l \2.4^ . In 
particular, the assumptions l \2.4\j are satisfied for following initial data 

Po(x) = (|x-xon^/^^°~^\ mo(x) = 0, xen. 

Without the loss of generality, the total initial mass is renormalized to he one throughout 
the present paper, i.e., 

Po{x)dx = 1. 



The boundary conditions are one of the boundary conditions of Dirichlet type and 
periodic type for Eqs. f l2.ip - fl2.2l) imposed as 

(1). Dirichlet case: 

pu{0, t) = pn(l, t) = 0, t>0, (2.6) 
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(2). periodic case: 

p, u are periodic in x of period 1, (2-7) 
where we consider Eqs. fl2.1l) - fl2.2p on M x (0, oo). 

Remark 2.2 For the case of Dirichlet boundary (\2.6[\ . the boundary is given by the 
physical observable momentum instead of velocity. It is natural to employ such boundary 
condition as considering the dynamics of (global in time) weak solutions to the IBVP 
for the compressible Navier-Stokes equations with possible vacuum states included since 
it is usually the case that, at vacuum states, the momentum is zero and is observed and 
controllable, but almost nothing is known yet for the velocity for weak solutions. Note 
here that for any (weak) solution away from vacuum at the boundary, the boundary 
condition {\2.6^ reduces to u{0,t) = u{l,t) =0, t >0. 

In order to define the weak solutions to the IBVP for the compressible Navier-Stokes 
equations f l2.ip - fl2.2l) with initial data fl2.3p and boundary condition (12. 6p or (12. 7p . we 
define the set of test functions as follows, 



and 



C^{n X [0,T)) for the Dirichlet case (1231) . 
C;^,(R X [0, T)) for the periodic case fl^Tp . 



C^{n X [0, T)) for the Dirichlet case fl^ . 
C~,(M X [0, T)) for the periodic case fl2Tp . 



with C;^,(M X [0,r)) defined by 

C^eM X [0,T)) = e C°^(M X [0,r))| ^ is periodic in x of period 1}. 

We define the weak solutions to the IBVP for the compressible Navier-Stokes Equa- 
tions dlllD-O as follows. 

Definition 2.3 (global weak solutions) For any T > 0, [p, u) is said to be a weak 
solution to Eqs. (\2.1^ - (\2.^ with initial data {\2.3^ and boundary value {\2.^ or {2.7) in 

nx{o,T), tf 



(2i 



< p G L-(0,T;Li(fi) nL^(l])), (p"-^/^)^ ^ L^{Q,T-L\n)), 
^pa G L-(0,T;L2(fi)), p^u^ G L\Q,T-W^^l'\n)), 

and (p, u) satisfies 

Poil>{x,0)dx + / / pt/jtdxdt + / / y/py/pwjjxdxdt = (2.9) 
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for any ^/^ G and 

f-T 



mo(p{x, 0)dx + / y/p{^/pu)(ptdxdt 
Jo Jn 

+ [ [ {{^/puy + p^) ^.Jxdt - {p^u^, if^) = (2.10) 



'0 Jn 

for all (fi ^ ^. The nonlinear diffusion term p°'Ux is defined as 

{p"u^,(p) = - j j p"'~^^'^y/puip^dxdt 
2a 



^ j {p''~^'\y/pu^dxdt (2.11) 



2a - 1 Jo 

for any </9 G where p G L°^{VL x (0,T)) due to ( \2.8^ . Moreover, for the spatial 
periodic case ( 2.7\) , (p, ^/pu) is also periodic. 



Remark 2.4 For the Dirichlet case, i\2.9^ . together with the fact p G L°°{Q x (0,T)) 
due to ( \2.8h . implies that {p,u) satisfies the Eq. I\2. 1} in the sense of distribution and 
justifies the boundary condition ( \2. gl) in the sense that, for any time interval I C [0, T], 

e~'^ / / \rP\fP'^{^i s)dxds — > 0, e^^ / / ^^uix., s)dxds 

Jl Jo Jl Jl-e 

as e ^ 0+. If further y/p^/pu G Lp(0,T; W^''J{n)) for some p > l,g > 1, then (El) 
yields that G L'p{0,T] Wq'''{Q)), that is, the Dirichlet boundary condition i \2. 6h 

is satisfied in the sense of trace. 

Definition 2.5 (global entropy weak solutions) Let {p,u) be a global weak solu- 
tion (in the sense of Definition \2. 3\} to ( \2.1\} - ([27^ with initial data l \2.3\} and boundary 
value ( \2.6i} or [2Sl) in Q x (0,T). Then, {p,u) is said to be a global entropy weak 
solution if there exists some function A G L^(f2 x (0,T)) satisfying ( \2.11\\ . i.e., 

J J Aifdxdt = — J J p°'^^^'^^/pulpxdxdt 



2a 



^ j {p''~'/^)^^/puipdxdt (2.12) 



2a - 1 

for any G $, and the following uniform entropy inequality holds 

sup / {\^u\^ + \ip''"'^^)^\'' + 7r{p))ix,t)dx 
o<t<T Jn 

+ f ! (i(p(^+"-')/').r+A')(x,t)M 

Jo Jn 
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<Co / + \{pr'%\' + Mpo)){x) dx (2.13) 

wt/i Co > independent ofT, and 

l:pTP^ ^/ 7>1, "^^^^ 1:^P^ ^/ 7 > 1- ^ ^ 
We have the following result on the existence of global entropy weak solutions. 



Theorem 2.1 (Global existence) Assume that 



1 

2' 



a 

7 > 77- 



(2.15) 



and that the initial data {po,mo) satisfies and E L^{il) for some positive 

Po 

constant v. Then for any T > 0, there exists a global entropy weak solution (p, u) to the 
IBVP for the compressible Navier-Stokes equations i\2.1\i - i\2.3^ with boundary condition 
or ( [^ inQ X (0,T) in the sense of DefinitionlKE 
Moreover, for the case of Dirichlet boundary condition i \2.6^ . if a G (1/2, 3/2) and v 
satisfies ^2.11^ {see Remark \KR below), then in addition to (l^.gp . the solution (p, ^/pu) 
satisfies 

^i^pu) G L2(0,T;<'(^+^'^)/(^+'^)(fi)), (2.16) 
i.e., ^/p{^/pu) satisfies the Dirichlet boundary i\2.6i} in the sense of trace. 



Remark 2.6 (1). Theorem \2.1\ above holds for the compressible shallow water equation 

(2). For the Dirichlet case {\2.6^ . one of the available ways (available for the case 
a G (1/2,3/2) and -1 > I) to obtam that ^^~pa G Lp(0, T; l^i''?((])), for some p > 
l,g > 1, is to choose the positive constant v in Theorem \2. 1\ such that 

/oraG(i,l], 7>1, 

/or a G(l, I), 7e[l,i±^), (3.17) 
for a e (1,1), 7>^- 



z/ G (0,27-a 
< u G 

' 4(a^l) 



2(a+7-2) 2(27-a) 
3— a— 7 ' 27 



G 



3-2a ' 



00 



Next, we show that there is a global entropy weak solution (p, u) in the sense of 
Definition 12.51 for which the vacuum states and the structure of interface, if existing 
initially, can be preserved for a short time, so long as the initial data has additional reg- 
ularity besides (12.41) and the fluids and the vacuum states in initial data are connected 
"smoothly" . In addition, the weak solution (p, u) is actually a unique regular solution 
for the short time. For simplicity, we consider the case of one point vacuum state 
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contained at x = Xq E (0, 1) in the initial data (po^^^o) = (PO)Po^o) with additional 
regularity 

Ao\x — xo\'^ < po{x) < Ai\x — xqI'^, for any x G fi, (2.18) 
uo e C\n), {pV^"'% e L'^iQ), Po'^'/'^(p>o.). e L''{n), j = l,n, (2.19) 

with n > 2 an integer; and in the case of continuous vacuum state of one piece initially 
on = [xo,a:i] C (0, 1) in the initial data, we require 

Ao{xo - xY < po{x) < Ai{xo - x)", X e [0,Xo), 
Po(x) = 0,mo(x) = poUoix) = 0, x e [xo,Xi], (2.20) 
Bo{x - xiY < po{x) < Bi{x - xif, X e (xi, 1] 



and 

(pr'^'^'')x e L'^in), J = i,n, Uo e c\n \ n': 



Po'^'^''{p>o.). G L'^i^ \ ^'), J = 1, 



(2.21) 



n, 



with n > 2 an integer. Here, a, Aq, Ai, and Bq, Bi are positive constants, and the 
power a G (a_,cr+) with positive constants a± given in (12.331) later. We also require 
that the initial data (po,mo) = (po)Po^o) given by (12.31) is consistent with boundary 
value for the Dirichlet boundary condition. 

We have the following results on short time structure of global entropy weak solu- 
tions. 

Theorem 2.2 (Short time structure of vacuum states) In addition to the as- 
sumptions of Theorem \2.1\ assume further that 

a > |, 7 > max{l, a) (2.22) 

and that there is either one point vacuum state in initial data {pq,Uq) with { \2.18^ - 
{ \2.iy^ satisfied or a piece of continuous vacuum states in initial data (po, ^o) 
{ \2.20^ - i\2.21\i satisfied. Then, there exists a global entropy weak solution {p,u) to the 
IBVP for the compressible Navier-Stokes equations ( \2.1^ - i^2^ with initial data ( \2. 3^ 
and boundary value (1^.61) or [2/1] in the sense of Definition \2. 5l 



Moreover, there is a short time T^ > 0, so that the global entropy weak solution 
{p,u) is uniqu^ and regular on the domain x [0,T,,], and the initial structure of 
vacuum states is maintained for t G [0, T*] in the following sense: 

For the case of one point vacuum state initially, {\2.18^ . the solution (p, u) is regular 
and unique on the domain Q x [0,T*], 

(p,m) G C°(0 X [0,r,]), G L°°(0,T,;C°(O)), (2.23) 



^Here the uniqueness is specified for density p and momentum pu = ^J^^u for continuous vacuum 
states of one piece. 
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ll^llL°°(nx[0,T.]) + ll^x||L°c([0,T*];CO(n)) < C(T^)- (2.24) 

The one point vacuum state propagates along particle path, namely, there is one particle 
path X = Xo{t) : [0, T,] ^ U with Xo{t) e C([0, T,]) defined by 

Mt)=uiXoit),t), Xo(0) =Xo G (0,1), (2.25) 

so that 

a_|x-Xo(t)r < p{x,t) < a+|x-Xo(t)r (2.26) 

for {x,t) G r2 X [0,T*], where the two positive constants a± are independent of time T^. 

In the case of a piece of continuous vacuum states initially, { \2.20^ . there are two 
particle pathes x = Xi(t) : [0,T*] Q with Xi{t) G C([0,T*]),i = 0, 1 defined by 

X,{t)=u{Xi{t),t), X,(0) =x, G (0,l),z = 0,l, (2.27) 

so that it holds for some positive constants a±, b± independent of the time T^,, 

a_(Xo(t) - xr < p{x, t) < a+(Xo(t) - x)^ (2.28) 

for {x,t) G [0,Xo(t)) X [0,T,], and 

b.{x-Xi{t)y < p{x,t) < b+{x~Xi{t)y (2.29) 

for {x,t) G {Xi{t),l] X [0,T*] respectively, and the interfaces separating the fluid and 
vacuum coincide with the particle pathes 

p{x,t) = 0, pu{x,t) = 0, {x,t) G [Xo{t),Xi{t)] X [0,T,]. (2.30) 

The solution (p, u) is regular and unique up to the vacuum boundary 

pec\nx[o,n]), M G c°(fi X [o,T*] (2.31) 

lkllL°°(nx[0,T.]\f75,J + ll'"2:|L°o(f7x[0,T,]\nO J ^ C{T*)- (2.32) 

where fi?,^ = (Xo(t), Xi(t)) x [0,T,]. 

Remark 2.7 (1). The constant exponents a± are defined as (j± = /3±/(l — (3±) > 
with (3± determined by 

(3. = max{^, i(l - I)}, = min{l, i(l - ^), ^(4 - 1)}, (2.33) 

while the positive constants a± are independent of the time T* . 

(2). The regularity assumptions (\2.1^ - (\2.19^ are satisfied for the following initial 
data 

Po(x) = |(Ao + Ai)(|x - a;o^)('^-+'^+)/^ uo{x) =0, xen, 
and the regularity assumptions (\2.2C^ - i\2.21\i are satisfied for the initial data 

i(Ao + Ai)(xo-x)('^-+-+)/2, XG [0,Xo), 
Po{x) = {0, xe{xo,xi), uo{x) = 0, xe^l. 

UBo + Bi){x - xi)('^-+'^+)/2, X G (xi, 1], 
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Next, we prove that for any global entropy weak solution (p, u) to the IBVP (12. ip - 
(12. 3p together with boundary condition (12. 6p or (12. 7p in the sense of Definition 12.51 
even though in some cases that the vacuum states may exist for some finite time, for 
instance, in the cases as shown by Theorem l2.2l any possible vacuum state has to vanish 
within finite time after which the density is always away from vacuum. Simultaneously, 
not only can the velocity field be defined in terms of the density and momentum, and 
the nonlinear diffusion is represented in terms of the density and velocity, but also the 
global entropy weak solution (p, u) is shown to be a unique and strong solution after 
the vanishing of vacuum states. We have the following result. 

Theorem 2.3 (Vanishing of vacuum states) Assume that 

a>^, 7 > 1. (2.34) 

Let {p,u) be any global entropy weak solution to the IBVP i \2.1\i -( ^27^ with initial data 
c/p and boundary value {\2.6} or 7p in the sense of Definition \2.5\ Then, there exist 



some time Tq > (depending on initial data) and a constant p_ so that 

inf p(x, t) > p_ > 0, t > To, (2.35) 

and the global entropy weak solution (p, u) becomes a unique strong solution (p, u) for 
t >Tq and satisfies 



p e L^{To,t; H\n)), pt e T-(To, t; L\9)), 
u E H\To, t- L\n)) n T2(To, t; H\n)), 

with velocity u and nonlinear diffusion term given by 



(2.36) 



n^^, (p"n,.)x = Ax, (2.37) 

respectively. In addition, for 

for the Dirichlet case, 

m^dx for the periodic case, 

there exist two positive constants po,Co both depending on initial data (po, mo) and p_, 
such that 

Up- Po,u-Us)i;t)\\L2^n)<Coe-'''^'^^'\ t > Tq, (2.38) 
where and what follows f denotes the average of f over the bounded domain Q, i.e.. 




f = ]^\ I fi^)dx = J f{x)dx. 
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Remark 2.8 (1). The Theorem \2. 3\ shows that any possible vacuum states must vanish 
in finite time. This theory applies to the compressible shallow water equation f I j . ^ . 
(2). It is easy to verify (see the proof of Proposition [377]) that the phenomena of 



vacuum vanishing (\2.35tj in finite time actually happens for any global weak solution 
{p,u) to the IBVP i \2.1h - l[K3i with boundary condition ( \2.6fi or i \2. 7p in the sense of 
Definition \2. 51 satisfying the following entropy inequality 

sup / (lv/pn|2+|(p"-V2)^|2 + vr(p))(x,t)dx+ r||(p(^+°-iV2)^j|2^^ 
o<t<TJn Jo 

<Co I + \{pr''\? + ^^{p.)){x)dx (2.39) 

with Co independent ofT. 

Finally, for any global entropy weak solution (p, u) to the IBVP (12. ip - (12. 31) together 
with boundary condition (12.61) or (12.71) in the sense of Definition 12. 5[ the density is 
continuous, i.e., p G CipL x [0,r]) for any T > 0, due to and ([23]). Thus, the 
continuity of p and Theorem 12.31 imply that if the density contains vacuum states at 
least at one point, then there exists some critical time Ti G [0, Tq) with Tq > given 
by (I2.35P and a nonempty subset fi*^ C such that 

'p(x,ri) = 0, VxGfi° 

p(x,Ti)>0, VxGf]\f]°, (2.40) 
^p{x, t) > 0, V (a;, t)eVtx {Ti,To]. 

It follows from (12.361) easily that for any S > 0, it holds 

-To 

\\ux\\L°°ds < oo. (2-41) 

Under the condition that vacuum states appear, we shall prove that the spatial deriva- 
tive of velocity (if regular enough and definable) blows up in finite time as the vacuum 
states vanish, even if the solution is regular enough for short time so that the velocity 
field and its derivatives are bounded as shown by Theorem 12.21 

Theorem 2.4 (Finite time blow-up) Let {p,u) be any global entropy weak solution, 
which contains vacuum states at least at one point for some finite time, to the IBVP 
for the compressible Navier-Stokes equations {\2.1\i - {\2.3^ with boundary condition { \2. 6^ 
or l\2. 7p in the sense of Definition \2. 51 Let Tq > and Ti G [0,To) be the time such 



that l \2.35\\ and ( 2.40^ holds respectively. 



Then, the solution (p, u) blows up as vacuum states vanish. Namely, for Ti satis- 
fying p.^Op and for given any fixed rj > 0, it holds 



lim / ||^^a||L°°c^'5 = oo. (2.42) 



t- 
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On the other hand, if there exists some T2 G (0, Tq) such that the weak solution (p, u) 
satisfies 

||w||Li(0,T2;VKi'°°(n)) < OO; 

then, there is a time T3 G [T2,To) so that the blowup phenomena happens for {p,u), 
i.e., 

lim / ||n2:||ioo(is = 00. (2.43) 
t-*T.- Jo 



Remark 2.9 Theorem 2.4 implies that for any global entropy weak solution (p, u) to 
the IBVP for the compressible Navier-Stokes equations { \2.1\i - {\2.2^ with initial data 
i \2. 3^ and boundary value ( \2. ^ or ( |^. 7p in the sense of Definition \2. 51 which contains 



vacuum states at least at one point initially, the finite time blowup phenomena i \2.42^ 
happens for such solution (p, u 



Remark 2.10 Theorems \2. l\^^ provide a complete dynamical description on the van- 



ishing of vacuum states and blow-up phenomena for the global entropy weak solutions 
to the compressible Navier-Stokes equations with density- dependent viscosity. That is, 
a global entropy weak solution exists for general large initial data with finite entropy. 
For short time, such weak solution is unique and regular with well-defined velocity field 
subject to additional initial regularity, and any existing vacuum state is maintained with 
the same interface structure as initial. Then, within finite time the vacuum states van- 
ish definitely and the velocity blows up {even if it is regular enough and definable along 
the interfaces). After the vanishing of vacuum states, the global entropy weak solution 
becomes a strong one and tends to the non-vacuum equilibrium state exponentially in 
time. This dynamical phenomena is quite similar to those well-known for the 3-D in- 
compressible Navier-Stokes equations. However, before the time of vacuum-vanishing, 
the uniqueness of the global entropy weak solution to the compressible Navier-Stokes 
equations with density- dependent viscosity subject to the initial data is not known yet. 



Remark 2.11 All theories established in Theorems{KW2^fit the shallow water equa 



tion (\1.2^ . We believe that such phenomena described by Theorems \2.1\ 2.4 are also 



observed for other compressible fluids with density- dependent viscosity, such as Navier- 
Stokes equations with capillarity and/or drag friction, Navier-Stokes-Poisson system, 
etc. 



Remark 2.12 It is interesting to investigate the (global) dynamics of (one- dimensional) 
interface connecting vacuum from initial time until the vanishing of vacuum in order to 
investigate the dynamics of interface and the vanishing of vacuum state and to verify 
the formation of singularity for general case. 

Remark 2.13 (1). Another interesting problem is whether the phenomena of the 
vanishing of vacuum states and blow-up happens for multi- dimensional compressible 
isentropic Navier-Stokes system with density- dependent viscosity, especially on spatial 
bounded domain with the Dirichlet boundary condition fl^.gp . 
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(2). It is also interesting to study whether any vacuum states shall vanish within fi- 
nite time and blow-up phenomena happens for the multi- dimensional full Navier-Stokes 
system with density- dependent viscosity. It is not obvious yet [although we expect) since 
it is not clearly understood yet how the dynamics of temperature and heat- conduction 
shall affect the global existence of weak solutions and the evolution of vacuum states, 
especially in the case of spherical symmetry under the Dirichlet boundary condition 
i \2. 6^ . This is under further investigation [26]. 



3 Global existence of entropy weak solutions 

In this section, we will establish the existence of global entropy weak solutions to 
the IBVP for the compressible Navier-Stokes equations (12. ip - ( 12731) together with f l2.6p 
or (12. 7p . Since the compactness arguments are straightforward in the framework of 
Mellet-Vasseur [38], we only need to construct a sequence of approximate solutions by 
using some ideas developed by Jiang-Xin-Zhang in [23] for one-dimensional compress- 
ible Navier-Stokes equations with free boundary, to establish some uniform a priori 
estimates, in particular, the lower and upper bounds of the density for the approx- 
imate solution sequence with the help of BD entropy, and to justify the boundary 
condition (I2.16P for (limiting) weak solution for the Dirichlet follows. 

Outline of proof of Theorem \2.1\ Step 1. Construction of smooth approximate so- 
lutions. Let us consider the following approximate compressible Navier-Stokes equa- 
tions inspired by Jiang-Xin-Zhang in [23] 

pet + {peUs)x = 0, (3.1) 
{peUe)t + (peU'l + Pip6))x " iPeipe)Uex)x = 0, (3.2) 
{pe,peUe){x,0) = (poe, "Iqe) (x) , (3.3) 

with one of the following boundary conditions 

u,{0,t) = u,{l,t) = 0, (3.4) 

or 

Pe,Ue are periodic in x of period 1. (3.5) 
The viscosity p^ is given by 

Peip)=p'' + ep', £>0, (0,1/2), (3.6) 

where we remark (see Remark 13.11 below) that although the modified viscosities (13. 6p 
are limited to the case 6 G (0, 1/2), the methods adapted here to construct approximate 
solutions also can be applied to the case 6 G (0, 1) and ^ > 9 after some modification, 
the interesting reader can refer to [16,23] and references therein, we omit the details 
here. The initial data posi^oe G C°°{VL) satisfies 

poe Po in L\Q), po7^/^ p^"^''^ in H^{n), 
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as e —>■ 0"*", and 

Poe > cos'/^'-''^ (3.8) 

for some Cq independent of e. By (13.71) we can assume without the loss of generahty 
that 

Posdx = 1. (3.9) 



By the standard arguments (see [23] for reference), after applying the classical theory 
of parabolic and hyperbolic equations, one can obtain that there exists some T* > 
such that the approximate problem f l3.ip -( l373l) together with the boundary condition 
(13. 4p or (13.51) has a unique smooth solution {pe,Ue) on [0,T*] with the density away 
from vacuum, i.e., 

ps{x,t) > 0, for all (x,t) G H x [0,T,]. (3.10) 

Step 2. The a-priori estimates. To extend the local solution globally in time, one 
needs to control the lower and upper bounds of the density and get some a-priori 
estimates. As mentioned in the introduction, this relies on the BD entropy inequality 
developed by Bresch and Desjardins (see [2,3,6,7] for instance). In the Eulerian 
coordinates, the BD entropy for the approximate solutions (p^, u^) of the IBVP problem 
for the system (I3.ip - (I3.1I) reads (see [2,3,6,7] for instance) as 



1'^ 








1 1 





+ r / {Pe{Ps)\ne.? + '''^''^f^'''\ pe.\^)[x,s)dxds < C, t>0 (3.11) 
Jo Jn " 

with positive constant C independent of t and e. Indeed, multiplying (13.21) by u^, 
integrating by parts on [0, 1], and using (13. ip . we get the classical entropy estimate 

/ {\y/p'eUe\^ + 7r{p,)){x,t)dx+ [ f fie{pe)\Usx\^ dxds < C , t>0 (3.12) 

Jn Jo Jn 

with positive constant C independent of t and e. 

Multiplying the transport equation (13. ip by (f'^ipe) with pf'ir{p) = Pe{p), i-e., ^e{p) = 
aP°' + as in [8], we get 

^eipe)t + Ve{pe)xUe + (p'.ipe) PeUex = 0. (3.13) 

Differentiating above equation with respect to x and denoting = ipe{.Pe)x/ Pe^ we have 

{.PeVe)t + {.PeVeUe)x + {.Pe{.Pe)Uex)x = 0. (3.14) 

Summing the equations (I3.14p and (13.11) together, multiplying by Ue+Ve and integrating 
by parts on [0, 1], we have after a straightforward computation that 

{\y^eUe\' + \^\' + 7lipMx,t)dx 

+ f f {pe{Pe)\Uex? + '-^^^^^\Pex?){x, s) dxds <C,t>0 (3.15) 

Jo Jn 
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with positive constant C independent of t and e. This together with (13.121) 

<2(pf-^ + ^Vf-^)(p,.)^ 

gives rise to the BD entropy (13.111) . 

In order to derive the lower bound and upper bounds of density for the approximate 
solution {pe,Ue)j it is convenient to make use of the Lagrangian coordinates {y,t) with 

Ps{z,t)dz,t = t, 

because the approximate solution (pcUs) is away from vacuum at least for short time 
due to (I3.10p and fluid density transports along particle path. In Lagrangian coordi- 
nates, we have the equivalent system for {ps,Ue) as follows 

Pet + pluey = 0, (3.16) 

Uet + ipipe))y " ipePeipe)Uey)y = 0, (3.17) 
{pe,Ue){y,0) = {poe,moePoe)iy), (3.18) 

together with the boundary conditions similar to ( 13. 4^ or ( 13. 5p . where (y, t) G x 
(0,T*] and 



fiz, = (0,Le), with L^= I poe{x)dx = 1 

Jn 



due to (Km . 

Correspondingly, the BD entropy equality ( 13. lip becomes in Lagrangian coordinates 
as follows 



sup E^{t)+ / / pePe{pe){Ueyfdydt 
0<i<T* Jo Jul 

+ r I ('^^'""^' + ^Z^^'"') (Peyfdydi < Coi^.(O) (3.19) 
with Cq independent of and e > 0, and Ei,{t) defined as 

Ee{t) ^ [ {u! + M)yf + e\{pl)yf + p-Mps)) dy. 

Since 

I p~^dy = I dx = \Q\ = 1, (3.20) 

it follows from the continuity of pe that there exists some ?/o(^) ^ such that 

p,{yo{t),t) = \QL\im = l. 
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Hence, it implies from ( I3.19P that 

p^(2/,t) = p^(2/o(t),t)+ r {p^)ydy 

Jyo{t) 

< {\^lLm~^T + CoE,{S)) + \^lL\. (3.21) 

This yields the uniform upper bound (w.r.t. time t G [0,T*]) for the density. 

To obtain the lower bound for the density p^, we employ the idea in [23]. That is, 
we consider the upper bound for = p~^, which can be estimated by (13.191) and (I3.20p 
as follow: 



Vs{y,t) < / Ve{y,t)dy+ / {Ve) \pey\dy 

<1 + C max(i;,)^+i/2 1| {p% ( / v,{y, t)dy) 

<1 + ^ maxt;, + C{e-^E,{0)y^^^-^'\ (3.22) 

This shows that for all {y,t) & x [0,T^], 

52/(1-20) 

Peiy,t) > 2^2/(1-20) +(^(^^(0))l/(l-2e)- (3-23) 

Thus, we have obtained the uniform upper bound (13.211) and the lower bound (I3.23P 
(w.r.t. time t G [0,T^,]) for density p^. Using these uniform bounds, one can extend 
the local smooth solution globally in time by standard arguments (see [23] for details). 
Moreover, the total mass of the approximate solutions is conserved due to the boundary 
conditions (13. 4p 

Pe{x,t)dx = I poe{x)dx, W t > 0. (3.24) 



Remark 3.1 Although the estimates i\3. 2^ and ( \3. 23^ are established for the approx- 
imate solutions in the case that the modified viscosities ( \3.(j^ are limited to the case 
9 G (0, 1/2), the methods adapted here can he applied to establish the lower bounds like 
( \3. 23^ in the case 9 G (0, 1) and ^ > 9. Indeed, one can also establish the following 
estimates, similar to [16,23], for the approximate solutions {p,Us) as 

sup /"(/)f(x,t)rfx<C(r) 

for any integer n > 1 and C(T) a constant independent of e,9. Thus, we can verify 
similar to i\3.22^ that 

v,<l + C [ vl+'\{p'),\dx < 1 + C(T)( / vf+^^dyy/" <C + C(T)tM^+i)-^)/'' 



H.-L. Li, J. Li, & Z. Xin 



19 



with q = 2n/{2n — 1). This together with Young's inequality implies the positive lower 

2n 



bound of density for integer n sufficient large hut fixed so that 9 < !- so long as 



(0,1). 

In addition to (13. lip , we show p^Ue is actually bounded in L°°(0, T; L^"*"'^(r2)) for 
some Tj > 1/(2 + u) with u > given by 

u > arbitrary if 27 — a > 1, 

i/e (0,2(27 -a)/(l + 0-27)] if 27-a G (0,1). 

Note here that for a G (1/2,3/2) the constant u defined by fl2.17p actually satisfies 
( K2^ . It follows from fl33T]) that 

sup II (pr 1/^)^1^2(0) <C, (3.26) 
t>o 

which together with fl3.24p and the mean value theorem gives rise to 

sup WpsWT-'^' < We''^'" + \\{pT''WWi.^) < C (3.27) 

where C denotes some generic positive constant depending on -E'e(O) but independent 
of T. With the help of f l3.27p . we can make use of the idea due to [38] to derive the 
following a-priori estimate 

sup / pe\ue\'^^''dx < C, (3.28) 

0<t<T Jo. 

which together with fl3.27p implies for 77(2 + z/) > 1 that 



/ {pl\u,\f^'' dx = [ pf+''^\u,\^+''dx<C sup / 
Jn Jq Q<t<T Jn 

and, in particular, for rj = 1/2 that 



p,|n,|2+^dx < C, (3.29) 



sup / iVP^Klf^'" dx <C (3.30) 
o<t<TJn 

with C one generic positive constant depending on both T and E^{0). 

Thus, we can follow the compactness arguments in [38] to prove the convergence 
of the approximate solution {p^^Ue) to some expected weak solution (p, m) as e ^ 0+. 
Namely, it holds 

Pe^P inC([0,T] xH), (3.31) 
{pT"^). - (p"^'/')x weakly in L\n x (0,T)), (3.32) 
^eUe Vpu, ^ p'^u in L^+^/^^l] X (0, T)), (3.33) 
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and there exists some function A G x (0,T)) such that 

py^x A weakly in L'^{Q x (0, T)), as e ^ 0. (3.34) 

Step 3. Justification of the Dirichlet boundary condition 1^2. 16\) for a G (1/2,3/2). 
For a G (1/2, 1], it follows from fl3TTD and (K29h that 

||(PeMe)a;.||i2(o,T;L(4+2<')/(4+^)(n)) 

<rlL°/2,, 11 I r ||('^"-V2^ ^l/2„, II 

^ HP's "e^-|lL2(0,T;L2{n)) ^ ^ II /x^e 1 1 L2(0,T;L(4+2'^)/{4+'-) (H)) 

< r-K^ ll/'n^-l/^N II . 11^1/2,, II 

- ^ II V/^e Jx\\L^{0,T;L2(n)) ll^e ^^W L^{0,T;L^+^{n)) 

< C. (3.35) 

For a G (1, 3/2), some additional estimates are needed here. First, for 7 > (1+q;)/2, 
it is easy to check that for z/ defined by f l2.17p that (3/2 — a)/ (2 + z/) > 1. Hence, for 
/? = (4 + 2i^)/(4 + v) one deduces from fl3:29l) and fl^:^ that 

\\{PeUe) x\\ L'^(Q^T;LP(Q.)) 

^ O ll/;^ "ez||^2(o,T;L2(Q)) ^ "-^ IKP'e j^Pe llL2(0,T;L/3(n)) 
-^^^\\\Pe /'a;llL°°(0,T;L2(f7)) ll^e ""^ 1 1 L°° (0,T;L2+''(f7)) 

< C. (3.36) 

In the case that 1 < 7 < (1 + a)/2, (12X71) yields that (3 - a - 7) (2 + z/) > 2. Hence, 
one derives from (13.291) and (13.111) that 

<rlL"/2,, II _L r ||('^{7+a-l)/2^ ^{3-a-7)/2„, II 

^ O ||/;^ "ex||^2(o,r.^2(Q)) ^ '-^ IH.P'e '^llL2(0,T;L'5{n)) 

< r -I- r II II lln(3-"-7)/2„, II 
^y^Ty^Wyp^ ;a;llL2(0,T;L2{C)) Ik£ 1 1 (0,T;L2+>'(f7)) 

< C. (3.37) 



Thus, the estimates ([MD, (l3:36|) and (13:3711 imply that (I2:T6|) still holds for a G (1, 3/2). 
The proof of the Theorem 12.11 on the existence of global entropy weak solution (p, u) 
is completed. □ 



4 Dynamics of vacuum states for short time 

4.1 Short time structure of vacuum states 

We prove the Theorem 12.21 in this subsection in order to study the short time 
structure of vacuum states for global entropy weak solutions. To this end, it is sufficient 
to show that there is a unique entropy weak solution in short time for the compressible 
Navier-Stokes equations ( 12.ip -( l2:2l) with initial data ( 12. 3p and boundary value ( 12. 6p or 
( 12. 7p under the assumptions of Theorem 12.21 as follows. 
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Proposition 4.1 (Vacuum states for short time) Under the assumptions of The- 
orem \2.2\f there is a short time > so that the unique entropy weak solution (p, u) 
in the sense of Definition \2.5i of the IB VP problem for the compressible Navier-Stokes 
equations l \2.1\} -l \2.2X} with initial data l \2.3\} and boundary condition l \2.6l} or l \2. 7] ) exists 



on the domain f2 x [0, T^] . The initial vacuum state (\2.19i) or ( \2.21t} is also propagated 
for the short time t G [0,T^'], more precisely, the properties {\2.23^ - {\2.26\} or l \2.21\j - 
i\2.32^ hold respectively for (p, u) . In addition, it holds that 

||(p"-^/').||L^(O.T^;L^(f.)) + IKWVp, {~puy+^/p'+n\\L^^O,Ti;LHn)) < C^'). (4.1) 

Proof: The proof of Proposition 14. II wiU be completed in subsections I4.2l - [4!3] later. 

Proof of Theorem \2.2i This is a consequence of the Proposition 14.11 and Theo- 
rem 12.11 In fact, the Proposition 14.11 shows not only that there exists a unique en- 
tropy weak solution (p, u) on the domain Q x (0, T^) in the sense of Definition 12.51 to 
the IBVP problem for the compressible Navier-Stokes equations fl2.ip - fl2.2p with initial 
data (12.31) and boundary condition (12.61) or (12.71) . but also that this short time entropy 
weak solution satisfies all the properties (I2.23I) - (I2.32I) . Now, choose time T^ = T^ — 6 
with (5 > a constant small enough. One can verify that at time t = T^, the density 
and momentum also satisfies the assumptions (12.41) . and particularly 

p(a;, T^) > on fi, p£t(.,T^) = 0, on {x E \ poix) = 0} , 
p{.,T,)eL'{n), {p--^i\.,T,%eL\n), (4.2) 

Thus, it follows from the Theorem 12.11 that there is a global entropy weak solution (p, u) 
for time t > T^, in the sense of Definition 12. 5l to the IBVP problem for the compressible 
Navier-Stokes equations (I2.ip - (I2.2I) with initial data 

(p,Pm)(x,T*) = {p,pu){x,T^), X eQ (4.3) 

and boundary condition (12. 6p or (12.70 . Define (p, u) as 

^^^N /(P'^)' {x,t)eQx[0,%], 

^ ' \ip,u), for (x,t) G X [0,r,]. ^ ■ ' 

It is easy to verify that (p, u) is a global entropy weak solution in the sense of Defini- 
tion [23] to the compressible Navier-Stokes equations (12. 10 - 02.20 with initial data (12. 3p 
and boundary condition (12. 6p or (12.70 under the assumptions of Theorem 12.21 This 
global entropy weak solution is actually unique and regular for short time t G [0, T*] 
and satisfies all the properties (12.230 - 02.320 . The proof of Theorem 12.21 is completed. 

□ 
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4.2 Compressible Navier-Stokes in Lagrangian coordinates 



In order to prove the Proposition 14. ![ we present an equivalent proposition for the 
compressible Navier-Stokes equations in the Lagrangian coordinates {y,t), instead of 
the Eulerian coordinates {x,t), through the coordinate transformation 



y = [ p{z,t)dz, xe{0,l), t>0 (4.5) 
Jo 

for the Dirichlet boundary condition, or 

l-Xo{t)+x 

y= p{z,t)dz, a;G(0,l), t>0 (4.6) 

Jxo{t) 

for the periodic boundary condition where x = XQ{t) is a particle path. In addi- 
tion, both the case of initial one point vacuum state and the case of a piece of initial 
continuous vacuum states will be studied with some additional initial regularities. 

We first describe the equivalent proposition for the IB VP problem for the compress- 
ible Navier-Stokes equations in the Lagrangian coordinates in the case of initial one 
point vacuum state. Thus, consider the compressible Navier-Stokes equations 

Pt + P^Uy = 0, (4.7) 
+ p{p)y - {pfJ'{p)uy)y = 0, y eQ, t> (4.8) 

with initial data 

(p, u)iy, 0) = ipoiy),Uo{y)), yen, (4.9) 
and the Dirichlet boundary condition 

u{0,t) =u{l,t) = 0, t>0 (4.10) 

or the periodic boundary condition 

(p, m) is periodic w.r.t. x of period one. (4-11) 

For the case of Dirichlet boundary, the initial data {po, Uq) is assumed to be consistent 
with the boundary values. Moreover, we assume that there is one point vacuum state 
aX y = yo for some fixed point yQ G (0, 1) and that the initial data is of additional 
regularity 

(4.12) 

uo G c\n), {pi^y, (p^"no,), G L\n) n L^^{n) 

with an integer n > 2, where Aq, Ai are positive constant, and the constant (3 G 
(/?_,/?+) with I3± determined by Remark 14. 3 [ Then the following result holds in the 
case of initial point vacuum state. 
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Proposition 4.2 (Point vacuum state for short time) Assume that (\2.2^ and 
i \4-lS^ hold. Then, there is a time > so that the unique regular weak solution (p, u) 
with point vacuum of the IB VP ( |^. 7p -( [7^ with boundary condition l \4-10^ or 1 14-1^ 

exists on the domain Q x [0, T^] and satisfies 

pec'{nx[o,T:])nc\[o,T:]-L\n)), (4.i3) 
u e c\n X [o,r:]) n c\[oX];L\n)), (4.i4) 

p^+X e L^{n X [0,T:])nC'/\[0,T:];L\n)), (4.15) 

||(p")2/IU°°([0,Ti];L2(n)) + \\pUy\\L^(nx[0,Ti]) < C X) . (4.16) 

In addition, the initial point vacuum state is maintained for the short time 

p(z/o,t) = 0, te[0,T:], (4.17) 

a-\y-yo? <p{y.t)<ci+\y-yo\^ {y,t) e [o,i] x [o,t,]. (4.18) 

Here, a± are positive constants independent of . 

Proof: The proof of the Proposition 14.11 will be given in subsection 14.31 

Remark 4.3 The choice of i3± > is such that 

p. = max{^, i(l - ^)}, = min{l, 1(1 - ^), ^(4 - i)} (4.19) 

for integer n > 2. It should be emphasized here that all the assumptions required here 
are satisfied for the shallow water equations U.S\) . i.e., 7 = 2, a = 1. 

Nest, we describe the equivalent proposition for the IB VP problem for the com- 
pressible Navier-Stokes equations in the Lagrangian coordinates for the case of initial 
continuous vacuum of one piece. In this case, we consider the IBVP problem for the 
compressible Navier-Stokes equation fl4.7p - fl4.9p with one of following boundary condi- 
tions: 



[1). mixed boundary condition 



M(0,t) = p(l,t) = 0, t>0, 

Ao{i - yf < Po{y) < A(i - yf, ye^, 



(4.20) 



or 

(2). mixed free boundary condition 



p(0,t) = M(l,t) = 0, t>0, 
5o/ < Po(y) < yen, 



(4.21) 



or 
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(3). free boundary condition 

p(0,t) = p(l,t) = 0, t>0, 

where, P G and Aq, Ai, Bq, Bi, and A'q, A[ are some given positive constants. 

The initial data is also assumed to be regular 

uo G c\n), {p2)y, (pr^Mo,), G L\n) n L^^{n) (4.23) 

with an integer n >2. Then the following results for free boundary value problems for 
the compressible Navier-Stokes equations hold. 

Proposition 4.4 (Continuous vacuum state for short time) Assume that (\2.2^ 
and { 4-23\) hold. Then, there is a time > so that the unique regular weak solution 



{p,u) of the IBVP problem for the compressible Navier-Stokes equations ( |^. 7| )-( [7^ 



with initial data ( |^.g| ) and free boundary condition ' \4-20^ , or ( |^.^jp , or l \4-22( ) exists 
on the domain Q x [0,T^'] and satisfies 

pec'{nx[oX])nc\[oX];L'm, 
ueC\Qx[0X])r^C\[0X];L'm, 
p'+"uy G L^{n X [oX])nc'/\[oX];L'm, 

and 

\\{P°')v\\L^{[0,Ti];L^{n)) + \\pUy\\L^{nx[0,Ti]) < CX)- 

In addition, the initial vacuum state is also maintained for short time, namely, it holds 
that 

a- (1 - y^ < piy, t) < a+(l - y^, {y, t) G (0, 1) x [0, T,'] 
corresponding to the mixed free boundary conditions f \4.2l^ , or 

b-y^ < p{y, t) < b+y^ (y, t) g (0, 1) x [0, T^] 

corresponding to the mixed free boundary conditions ' \4-21} , or 

c_(y(l - y)f < p{y, t) < c+(y(l - y)f, {y, t) G (0, 1) x [0, T,'] 



corresponding to the free boundary l \4.22!( ). Here, a±, b± and c± are positive constants 
independent ofT^. 

Proof: The short time existence, uniqueness and regularity of weak solutions 
of the free boundary problems for the compressible Navier-Stokes equations are well- 
investigated by many authors (see [23,32,33,56,57]). Proposition 14.41 can be proved in 
a similar way as [23,32,33,56], so we omit the details. 
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Remark 4.5 Based on Proposition \4-4\ foi" the compressible Navier-Stokes equations 
'^'^th either mixed free boundary condition l \4-20^ and l \4-21^ or free boundary 



l \4-22^ and the coordinates transformation from the Lagrangian coordinates to the Eule- 
rian coordinates, one can study easily the IB VP problem for compressible Navier-Stokes 
equations ( \2.1\) - ([27^ in the Eulerian coordinates with either the Dirichlet boundary con- 
dition h2.6\) or the periodic boundary condition < \2. 7p in the case of initial continuous 
vacuum state of one piece {\2.20i} . 

In fact, for the compressible Navier-Stokes equations l\2.1\i - l\2.2t} with Dirichlet 
boundary condition i\2.6[) in the Eulerian coordinates, the short time existence of unique 
solution subject to the case of a piece of initial continuous vacuum state {\2.2(Ji} - {\2.21\i 
in initial data can be constructed, in the Lagrangian coordinates, by combining two 
mixed free boundary value problems i \4-20^ and l \4-21\ ) together with one continuous 



vacuum state in-between as follows. Denote two particle pathes x = Xi(t) (assumed 
to be definable for short time ) starting from the initial vacuum boundary x = Xq and 
X = Xi respectively as 

X,{t) = u{X,{t),t), Xi{0) = x,, 1 = 0, 1, 
along which the vacuum boundary moves in the Eulerian coordinates so that 
{p,pu){x,t)=0, xe[Xo{t),Xi{t)], t>0, 

(4.24) 

pu{x,t)>0, X G [0,Xo(t)) U 1], t>0. 

We first choose the coordinate transformation from the Eulerian coordinates to the 
Lagrangian coordinates as 

y = Iq ^)'^^' ^ ^ [0' ^o(^)], 

X it) . ^^-25) 

yo = Jq " p{z,t)dz = Jq° po{z)dz < 1, conservation of mass, 



which gives y G [0,yo] and the mixed free boundary conditions { 4 -20^ . The application 



of the Proposition \4-4 with Q replaced by (0, yo) implies, via the inverse coordinate 



transformation x = p~^{z,t)dz, y G [0,?/o]; the existence of unique solution {pi,ui) 
of the compressible Navier-Stokes equations {\2.1\i - {\2.2i} on [0,Xo(t)] x [0,T=k] in the 
Eulerian coordinates with the initial data 

{pi,ui){x,0) = {po,uo)ix), xg(0,xo), po(a;o) = 0, 

and mixed free boundary conditions 

u{0,t) = 0, p(Xo(t),t) = 0, tG[0,r,]. 

Next, we choose the coordinate transformation from the Eulerian coordinates to the 
Lagrangian coordinates as 

y=l- J^p{z,t)dz, xG[Xi(t),l], 

(4.26) 

yi = j^^f^^^ p{z,t)dz = J^^po{z)dz < 1 conservation of mass, 
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which gives y E [1 — yi, 1] and the mixed free boundary conditions l \4.21\ ). Then Propo- 
sition implies, via the inverse coordinate transformation x = 1 — p^^{z,t)dz, 
y E [1 — yi, 1], the existence of unique {pr,Ur) of the compressible Navier-Stokes equa- 
tions l \2.1\\ - ([27^ on [Xi(t),l] X [0,T*] in the Eulerian coordinates with initial data 

{pr,Ur){x,0) = {po,Uo){x), xe{xi,l), po{xi) = 0, 

and mixed free boundary conditions 

M(i,t) = o, p(Xi(t),t) = 0, te[o,r,]. 

Consequently , we can construct the short time unique solution {p,u) to the IBVP 
problem for the compressible Navier-Stokes equations ^2.1\j - ^27^ with the Dirichlet 
boundary condition ( \2.(^ and a piece of continuous vacuum state (\2.20{) in the initial 
data as 

({pi,piui), on [0,Xo(t)] X [0,T,], 
{p,pu) = I (0,0), on (Xo(t),Xi(t)) X [0,T,], 

[{Pr,prUr), On [Xi (t) , 1] X [0, T^] . 

Similarly, we can obtain the short time existence of unique solution to the IBVP 
problem for the compressible Navier-Stokes equations l\2.1\} - l\2.2X\ with periodic boundary 
condition d^. 7| ) which can be viewed as a free boundary problem after the choice of the 
spatial reference point. The details will be omitted. 



4.3 Proof of Propositions I4.1H4.2 



We first prove the Proposition 14.21 in this subsection with the help of the a-priori 
estimates for (regularized) solutions and the construction of approximate solutions by 
a finite difference scheme, due to the modification of the ideas used in [23,33,43,56]. 
Without the loss of generality, we only prove Proposition l4.2l in the case of the Dirichlet 
boundary conditions with one point vacuum state in the initial data. 

First, we can easily derive some identities for (regularized) solutions as in [23,56]. 

Lemma 4.6 Let T > and assume that the solution (p, u) of the IBVP problem ( |/^. ?[ )- 
d^. iOD exists for t G [0,T] with p{yo,t) = 0. Then, under the assumptions of Proposi- 
tion \4.2\ it holds that 

ry ryo 
p^+'^Uyiy, t) = / ut{z, t)dz + p\y, t) = - Ut{z, t)dz + p\y, t), (4.27) 
Jyo Jy 

/"* /"* ryo 

p°'{y,t) + al p'^{y,s)ds =pQ{y) + a I j ut{z,s)dzds 

(4.28) 



ry 

--pQ{y)-a / / ut{z,s)dzds. 

Jyo 



for all y E Q with y ^ yo. 
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We also have the following useful a-priori estimates, whose proofs are similar to 
those for (13.111) and fl3.19p and thus will be omitted. 



Lemma 4.7 LetT > and assume that the solution {p,u) of the IB VP problem l\4-7\) 



( 1^ . j op exists for t G [0,T]. Then, under the assumptions of Proposition it holds 
that 

+ \\pmT\n) + fj \{p^'^''^'\?dxds < Co, 1 6 [o,r] 

||p(t)||L-(n) + ||(p°(t)),||L^(f^) < Co, t e [0,T] (4.29) 
with Co > a constant. 

Lemma 4.8 Let T > and assume that the solution (p, u) of the IB VP problem ( |/^. ?[ )- 
( \4-10^ exists for t G [0,T] with p{yQ,t) = 0. Then, under the assumptions of Proposi- 
tion \4-S\ there is a time T^ G (0,T] {depending on initial data) so that it holds that 

j Ut\y,t)dy + j u/~^P^^"ultdxds < Ci (4.30) 

uniformly for t G [0, T^'] with j = l,n and Ci > a positive constant, and that 
pUy{t) G C^^Cl) is uniformly bounded for any t G [0, T^'] 

\\pUy\\L^(^y^[i),Ti\) < C2 (4.31) 

with C2 > a constant. Moreover, the solution (p, u) is continuous 

(p,u)gC°(Ox[0,T:]), (4.32) 

and the initial one point vacuum state is maintained 

a- \y -yof< Piy, t) <a^\y- yo\^ (y, t) E n x [0, T,'] (4.33) 

with a+ > a_ > two constants independent of time T^ . 

Proof: Let us first assume that the weak solution is regular enough so that we 
can differentiate it through the equations and the interface as in [56], and the density 
is of the form 



a 



* b - l/ol^ < p{y, t) < a* |y - yol^ (y, t)e^x [0, T] (4.34) 



with a*, a* two positive constants to be determined later. It will be assumed further 
that the following a-priori estimate holds 

\\pUy\\L^(Qx[Q,T]) < Mo (4.35) 
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for some positive Mq to be determined later. We will prove (14.301) only for j = n below 
since the other case can be treated similarly. Taking inner product between (14.81) ^ and 
2n(ut)^"~^ over Q leads to 



d 
~dt 



ut\'^dy + 2n {p^)yt{utY^~^dy-2n {p'+'^Uy)yt{uty^-'dy = 



from which, one deduces after integration by parts and using Eq. fl4.7|] that 



d 



\ut\''''dy + 2n{2n-l) / p'+'^ul.iutY^-^dy 



=2n(2n-l)7 J p^-^ ptUytiutf'^-^dy 
+ 2n(2n - 1)(1 + a) / p'' ptUyUytiutf"-^ dy 



+ C I p'^-ul{mr~'dy 



(4.36) 



where one has used (12.220 . (I4.29p . the a- priori assumption (I4.35p . and Holder's inequal- 
ity. The last term on the right hand side of (14.361) can be estimated as follows. For the 
case a > 1, it follows from (14.291) and (I4.35P that 



<CMUI + \\ut 



2n+2\ 
L2" ) 



For the case a G (|, 1), since it holds by (14.270 that 



p'+XI < I / Ut{z,t)dz\+p^ < \\Ut\\L2. ■ |l/-yo|^'"-'^/'" + P^ 



(4.37) 



(4.38) 



which together (I4.34p implies 



(3+o)n 4n| 



IP-— =|p-(^+3°)"(pi+"nj,)4'^| =p 



-(l+3a)n 



Ut{z, t)dz + p'^ 



yo 



An 



<Cp-^'-^'-)-{\\ut\\% ■ \y - yoP^'"^'^ + p'"^) 

<C\\ut\\% ■ \y - y^|-(l+3a)n/3+2(2n-l) ^ ^3nh-a)+n(^-l) ^ 



(4.39) 



the last term in (I4.36P is estimated by 



p'^'^u'yiu,r--'dy\ < 



1/n 
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<C\\Utfj^,-'' ■ ( \\ut\\% I |y-yo|-('+'")"^+'('"-')rf2/+ / p3n(7-a)+n(7-l)^^ ■ 
\2n 



<C(ll«t|li^^' + l)- (4.40) 



where one has used the fact /3 < due to (ICTD . Substituting (14:371) and (ICTj) 

into fl4.36p shows 

1 1 \u,\'-dy + n{2n - 1) J p'-^'^ul.iu.Y^-'dy 
<CM',{1 + Mo^)(l + \\u,{t)rL^t%) < C{1 + M^){\\u,{t)rL^t' + 1). (4.41) 

Set 

T — ' r 2"-l ll"*(0)ll^2n //I An\ 

la — min|^2nC'(l+A/4)Ci ' C(l+M4) ' ^ J • I'^-'^'^J 

One can apply the Gronwall's Lemma to obtain (14.301) for t G [0, Ta] with Ci given by 

Ci =: 2||ui(0)||i^„ > ||n<(0)||i^„ + CT,(1 + Mq^). (4.43) 

To prove (I4.3ip and ensure the a-priori assumption (I4.35p . we use the equahty (I4.27P 
to get that near y = yo, it holds that 

puy{y,t) = p-''{y,t) r ut{z,t)dz + p^-'^{y,t), y ^ y^, (4.44) 
which, together with ([222]), ^AJM . and the fact /3 < ^ (l - ^) due to fjllB . implies 



2/ 

J.. ^ 



<Co'-" + - yor-'/'"-"''lk*(t)IU- 

<Co'-° + a^C^'" < CJ-" + 2ar||n,(0)|U2„ 
= : C2 < Mo 

so long as the constant Mq is chosen as 

Mo = 1 + CJ^" + 2aT\\um\\L^- (4.45) 
Next, we verify the a-priori assumption (I4.34p . Set 

T5 = min{^-^, TJ. (4.46) 

It follows from the equation (14.280 for y ^ yo and t G [0, T^] that 

p"{y,t) + a p"'{y,s)ds = pQ{y) + a I / ut{z,s)dzds 



"''2/0 



>^^o b - + \y- yor^(^A^ - Cl/'-aT,\y - y,\'~"'-'^^) (4-47) 
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where we have used the facts \y — yo\ < 1 and /5 < ^ (l — ^) due to fl4.19p . On the 
other hand, it follows from fl4.28p that 



rt rt i-y 

p°'{y,t) + a p^{y,s)ds = p^iy) + a 1 / Ut{z, s)dzds 

Jo Jo Jyo 

<A-,\y - yol"^ + at\y - yo\'-'^''lMt)\\L.. 
<K + Cy%t)|y-yor^. 
Define T, G (0,T,] by 



Set 



(4.48) 



Tc = mm{ — , T„T,}. (4.49) 



Zit)= [ p^{y,s)ds. 
Jo 

It follows from fHIlHD that 

{Z\t)r/'' <{AUCy'''^t)\y-yor^ 
which implies for t G [0, T^] that 

p^{y, s)ds <TMi + Cl'^''aT,r'^\y - 2/o|^^-"^^|y - yoV^ 
<CT,(A? + Cy'"aT,)^/"|2/ - 
As a consequence of (14.471) . (I4.50p . and (14.490 . one gets 

p"(y,t)> lA<^^\y - y,\-P - a f p\y.s)ds> \A^^\y - y,\-P . (4.51) 



(4.50) 



From (I4.48P and (I4.5ip . we can verify the a-priori assumption (I4.34p and justify the 
property (14.330 by simply choosing 



a* = a_ 



(A^/3)^/", a* = a+ = {A1 + C^^^aTj^/" (4.52) 



for any t G [0,T^'] with time T'^ = Tc determined by (I4.49p . One then derives from 
Eq. dSD, flCTj) . and fOTD that 

G L^{OX,H\n)), {p-)t G L^{QX.L\n)). (4.53) 

while flOTD . boundary condition flCTj) . and flCTj) for j = 1 imply that 

ueL^{OX,W^'''{m UteL'^iOX^L'm (4.54) 
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for any p G (1, /3 where one has used the fact P ^ > > I so that 

sup \\uy\\LP{n)= sup \\p'^\\mn)-\\pUy\\loo^^y,,oTn)<'^- (4-55) 

teio,n] teio,n] ' 

f l4.53p - p:.55p imply the continuity (14.321) of the solution {p,u), and the continuity of 
puy follows from the equation f l4.27p and that of (p, u) . The proof of the lemma is 
completed. □ 

Using Lemmas 14 . 6H4 . 8 1 and a direct computation, we can obtain the following result: 



Lemma 4.9 Let > be given in Lemma^T^ and (p, u) be the solution of the IBVP 
problem ( |/^ . 7[ ) -( [7^T0| ) . Then, under the assumptions of Proposition p^~^°'Uy{t) G 
0^(0,) is uniformly bounded for any t G [0,T^'] 

lim p^+"M,(t/,t) = 0, ||p'+"M,|Uoo(o,r-co(n)) < CiZ), 
y^yo 

and (p, u) satisfies for < s < t < that 

\\p{t) - p(s)|U2(^) + \\u{t) - u{s)\\L2^n) < C{T:)\t - si 
||p^+XW - p'-^'Uyis^L^n) < C{T:)\t - s\'/'. 

Proof: The facts that for any t G [0,T^], it holds that p^+'^Uyit) G C°{n) and 
limy^y^ p^~^°'Uy{y,t) = are due to the continuity of right hand side terms of fl4.44p 
and f l4.34p . By (14.3 ip and (I4.33p . one can check easily that 

\\P^^"Uy\\Loo^nxlO,Ti]) < \\p\\L°°{nxio,Ti]) • l|P^?/IU°°(nx[0,T^]) < C{Tl). 
Making use of Eq. (14.70 . (I4.30p and (14.310 . one can obtain 



and 



Wpi-t) - pis)\\L^n) < / \\pt{r)\\mn)dT = \\ / p Uy{T)\\L2^n)dT 

J s J s 

<C{t - s)||p||L°o(nx[o,r^]) ■ llp^lU-cnxicT^]) < C{Tl)\t - s\, 

\\u{t) - M(s)||i2(f^) < II / Ut{T)dT\\L2f^^) hth^iU) < C{Tl)\t - S\ 



p^+XW-P^X(^)IU^(f^) < II / {p'^''Uy{t))tdT\\L2^n) 



<C{\\ I Pip"%(r)|U.(n) + || / p'^^Uyt{r)dTU2^n)) 



<C{\\ / p^+XMIU^(^) + ll / P^+X*(r)rfrlU.(^)) 
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<C{T:)i\t - S\ \\p'+''ul\\L^^nx[0,Ti]) +\t- s\'/^p^'+^^/\t{T)\\L2^nx[0,Ti])) 

<C{T:)\t-s\'/'. 

The proof is completed. □ 



Proof of the Proposition 4-^ With the help of Lemmas I4.6H4.8[ we are ready to 



prove the Proposition | 

(1) Existence of weak solution for short time. We only deal with the case for the 
Dirichlet boundary condition and one point vacuum state in the initial data, the case 
of periodic boundary and one point vacuum state in initial data can be done in a 
similar way. Once the a-priori estimates are established as in Lemmas I4.6H4.9[ we are 
able to prove the Proposition 14. 2[ First of all, we construct a sequence of approximate 
solutions by modifying the finite difference scheme used in [33,40,43]. Without the 
loss of generality, we assume 

yo = ^, namely, po(|) = 0. (4.56) 

For any given positive integer iV = 2A; + 1 with A; > an integer, let h = 1/N. Consider 
the system of 2N ordinary differential equations 



d h f h \2^2n+2 ^ 

^P2n+1 + lP2n+l) ^ 



d H , (P2n+ir - (P2n-ir (4.57) 
dt""^- + h 

= ^{(pL+i)'+"«+2 - ul)/h - {pL.,Y^''{ul - Ul_,)/h}, 

with the boundary condition and point vacuum 

4 W = umit) = 0, p2fc+i(t) = 0, (4.58) 

and initial data 

P2nH-i(0)=Po((2n + l)|) <(0) = V(2n)|) (4.59) 
where n = 1,2,3, N. Here, we also assume 

P2iV+l(^) = P2^-lW' Plit) = P-lW' 4N+2it) = Ul,it) = 

which is consistent with the boundary condition and the fact that density is continuous 
and non-zero at the left boundary. 

By applying the idea as in [32, 43] and the similar arguments mentioned above, 
we can obtain the following uniform (w.r.t. h) a-priori estimates about the solutions 
(P2n+i5^2n) (here and below we omit the symbol h for simplicity) of fl4.57p - fl4.59p 
similar to Lemmas 14.6114.91 Details will be omitted (the reader can refer to [32,43] for 
similar arguments in details). 
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Lemma 4.10 Let {p2n+i,U2n) be the solution of i \4.5T\) -l \4-59^ . Then, it holds 

N t N 

E(l<W+^(/'2n+l(t)))/^+ / J2plnti r"'^V'"^'^ nds 
n=0 "^^ 71=1 

N 

= E(i<W + ^('^2„,+l(0)))/l, 

71=0 

N-l N-l 
71=0 n=0 

It follows from Lemma 14.101 and the standard theory of ordinary differential equa- 
tions that there exists a global solution (p2n+i 5^271) to fl4.57l) - fl4.59l) for any fixed pos- 
itive and h. Furthermore, the following properties hold: 

Lemma 4.11 Let (p2n+i; ^271) with n = 0, 2, 3, A^ — 1 with N = 2k + l be the solution 
of ( [737D -( |73PD . Then, zt holds forn>k + l that 



r d f 

P2n+iW =P2n+i(0) -a / ^ -f+^2j{s)hds - a / pln+i{s)ds, 

Jo +1 at Jo 

l+a (.^ U2n+2it) -U2n{t) _ ^ d 



j=k+l 

and for n < k that 



P2n-l{t) = P2n-M - a 



^ d /"* 

/ y^-rU2j{s)hds- a / pln-i{s)ds, 
Jo j^^at Jq 



Lemma 4.12 Under the assumptions of Proposition there is a short time T^ > 
so that it holds for n = 0, 1, 2, that 

N 

71=1 

a_|i(2n +l)h-\\^ < p2n(t) < a+\\{2n + l)h - \\^ , 
and for m = 1 or n that 

N „t Af 

J2{iu2,{t)r'h+ / Ep2;ri(^)(l^2.(3))^-^( ^"^"^^^T "^"-^^^V )/^ < c{t:). 

3=0 J^ j=0 

Here C{Tl) > and a± > are constants. 
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Lemma 4.13 Under the assumptions of Proposition and Lemma ^.l^ it holds for 
t G [0, Tl] that 

N 

\U2n 

n=l 

N 



l(P2n+l)'+"(n2n+2 - U^n) / h - (p2n-l)'+"(n2n " U2n-2) / h\ < C{Tl), 

n=l 

N N 

\P^nit) - P2n{s)\^h + Y \u2n^l{t) " U2n^3\^h < C(T:)|t - , 
71=1 n=l 

N 



n=l 



With the help of Lemmas l4.10l - [4!T3l we can define the sequence of approximate 
solutions {ph,Uh) on the domain Q x [0,Tl] as 

Phiy,t) = p2n+l{t), 

Uh{y, i) = \ [{y -in - \)h)u2n+i{t) + ((n + \)h - y)u2n-i{t)] 

for y e {{2n)\h,{2n + 2)\h). It can be verified that the following properties hold for 
the approximate solutions 

o / ,N U2n+l{t) - U2n-l{t) 
dyUh{y,t) = 

and 

\y-yo\'^ < phiy,t) < a+ly-yof, 

My,t)\ < C{t:), \p,{y,t)d,Uhiy,t)\ < C{Ti), 

for {y,t) G X [0,T^']. Then, the existence of weak solution for short time t G [0,t'^] 
follows from Helly's theorem, the diagonal process together with Lebesgue's theorem, 
and the a-priori estimates (see, for instance, [33]). The details are omitted. 

(2) Uniqueness of weak solution for short time. We prove the uniqueness of weak 
solutions for a G (|,7)- Without loss of generality, only the case for the Dirichlet 
boundary condition fl4.10p will be studied. Let (pi,ni) and {p2,U2) be two weak solu- 
tions of Eq. flI7n) - fl4.10|] satisfying Lemma SSHSl Denote 

n = pi-p2, iJ = Ui-U2, G (0, 1) X [0,T]. 

Obviously, the new unknown (n, ■?/;) with ip{0,t) = = n{yQ,t) = satisfies 



n 
P1P2 



+ = 0, yo, (4.60) 
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+ {Pj - Pl)y - (Pi+"^.). - ((Pi+" - P\^'')n2y)y = 0, (4.61) 

for {y,t) G (0, 1) X (0,T*] with zero initial data 

(n(|/,0),^(y,0)) = (0,0), ye (0,1). (4.62) 
Take inner product between p'^p'^^n and (14.601) over [0,?/o) U {.Voi 1] to obtain 

^ y Pi^^'^p2'^n^dy = -{l + a) j Uiyp'^p^'^n^dy -2 j p^p^^n^pydy 
<(C||pini,|Uo.+4) I p^'+y^'n'dy+]- f pl+^tfj^dy 



4 

1 
4 



<C(T,) ^ pr'+V2 + - ^ (4.63) 



where ^ fdx =: JJ"' fdx + J^^ fdx. Taking inner product between ip and ( 14.610 over f2, 
and noting that 'ip{0, t) = ip{l, t) = n{yo, t) = and pY^°'ui{yQ, t) = 0, i = 1, 2, we have 

^|llV'(t)f + jp\^''i^ldy = j{pl - pl)i,ydy - y"(p}+" - p\^^)u^y^l,ydy 

/p^ — p^ f p^"^" — p^^" 
-nil)ydy - (t — — nu2y^pydy 
Pl-P2 J Pi- P2 

<cj ' Pi'-n'dy + i j pl^'^r^dy 



< - / p{+>Jrfi/ + C pr^+>2 'ri^rfy. (4.64) 



(4.65) 



Summing up the two differential inequalities fl4.63p and (14.641) leads to 
f^mW + P^'^yfnit)dx + j>p\^^^ldy 

<C{T^) fp^'-'yfn'dy 
which, together with the initial data (I4.62p and the Gronwall's lemma, gives rise to 

WM^W + ^ Pi'^''P2^n\x,t)dx = 0, te [0,T]. (4.66) 

This together with pi{yo,t) = p2{yo,t) = implies the uniqueness of weak solution 

(pi,ui) = (p2,U2). (4.67) 
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The proof of Proposition 14.21 is completed. □ 



Proof of Proposition \4.1\ With the help of Propositions I4.2H4.4I and the inverse 



transformation from the Lagrangian coordinates to the Eulerian coordinates, we are 
able to prove the Proposition 14. 1[ Since the case of periodic boundary conditions with 
one point vacuum state in the initial data can be dealt with in a similar framework, we 
only deal with the IBVP problem for the compressible Navier-Stokes equations (12. ip - 
(12. 3p for the case of the Dirichlet boundary condition (12. 6p with one point vacuum in 
the regular initial data (12381) -( ICTI ). 

First, one can check easily that the IBVP problem (^Ji)-(^M, ([2l]), and fl238|) - 
(I2.19P in the Eulerian coordinates is equivalent to the corresponding IBVP problem in 
the Lagrangian coordinates for the equations (I4.7p - (l4.9p with the Dirichlet boundary 
conditions (14.10p and ( I4.12P through the coordinate transformation 



y = [ p{z,t)dz, xe(0, 1), t>0. (4.68) 
Jo 



Note that y G [0, 1] due to the conservation of mass. More importantly, the case 
of one point vacuum state in the initial data (I2.18p - (l2.19p is reformulated into the 
corresponding case ( 14.12P with yo = po{z)dz and /3 = It is easy to verify 

that all the assumptions of Proposition 14.21 are satisfied. Then Proposition 14.21 gives 
the existence and uniqueness of weak solution (p, u) satisfying (14. 130 - 04. 180 in the 
Lagrangian coordinate to the compressible Navier-Stokes equations (I4.7p - (l4.9p with 
the Dirichlet boundary conditions (I4.10p and one point vacuum state in the initial 
data. In particular, 

p,ueC\[0X];C'm, puy{t)eL°^ioX;C'm, (4.69) 

p{yo,t) = 0, p{y,t)>0, yj^yo. (4.70) 

These in turn imply, in terms of the inverse coordinate transformation of (14.680 . i.e., 

ry 

X = p-^{z,t)dz, y e [0, 1], 
Jo 

the existence and uniqueness of weak solution (p, u) on the domain (x, t) G [0, 1] x [0, T^] 
to the compressible Navier-Stokes equations (I2.ip - (l2.3p with the Dirichlet boundary 
conditions (12.60 and initial one point vacuum state (I2.18p - (l2.19p . Moreover, one can 
check that there exists one particle path x = Xo{t) defined by 

Xo{t) = u{Xo{t),t), t > 0, Xo(0) = xo 

satisfying 

yo= piz,t)dz = / poiz)dz 
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due to the fact 



d 



y{Xo{t),t) 



p{z,t)dz = 0, t>0. 



dr' ' dt 

This, together with (14.701) . gives rise to 

p(Xo(t),t) = 0, and p(x,t) > 0, x ^ Xo(t), t > 0. 

Moreover, it is easy to verify that the solution (p, u) satisfies all the properties (I2.23p - 
fl2.26p . and particularly 

p(x, T*) > on m{.,T^) = pu{.,T^) = 0, on {x eVL \ poix) = 0}, 

p(.,T,) e Li(fi) nL^(fi), (p"-V2(.,T,)), G L'in), 
e L\n). 



/5i+'^(x,T*) 



The case of periodic boundary can be treated in a similar way, we omit the details. 
The proof of Proposition 14.11 is completed. □ 



5 Vanishing of vacuum states and blow-up phenom- 
ena 

We shall prove that for any global entropy weak solution (p, u) to the IBVP for 
compressible Navier-Stokes equations (12. ip - (12. 31) together with (12.60 or (12.70 . any pos- 
sible vacuum state vanishes in finite time and the velocity (if definable and regular 
enough) blows up in finite time if vacuum state appears, for example, the density con- 
tains vacuum initially as in Theorem 2.2. The weak solution becomes a strong one 
after the vanishing of vacuum states and tends time-asymptotically to the non- vacuum 
equilibrium state exponentially. 



5.1 Vanishing of vacuum states in finite time 



Proposition 5.1 Let \2. 34\ ) hold. For any global entropy weak solution {p,u) to the 
IBVP for compressible Navier-Stokes equations (1^. ip -f l^^^ with initial data - {\2.3^ and 
boundary condition ( \2. 6^} or [2/1) in the sense of Definition \2.6l there exists a time 
Tq > such that 

inf p(x,t) > 0, for all t > Tq. (5.1) 



Proof: To prove (15.10 . we will employ an idea which has been used in [28] (see 
also [22, 27]) to show the blow up behavior of both the global strong solutions to the 
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IB VP for fl2.1l) - fl2.2p with the constant viscosity and the global strong solutions to the 
Stokes approximation equations, with initial data containing vacuum states. 

Let T G (0, oo) be fixed. In this subsection, C denotes some generic positive 
constant independent of T. First, it is noted that the total mass is conserved the total 
for any t G (0, T] 



p{x,t)dx = / po{x)dx 



(5.2) 



Based on the entropy inequality fl2.13p . it can be deduced from fl5.2p and (12.390 that 
for a constant b > max{a + 7 — 1, 2q; + 1, 1}, 



0<t<T 

1 1 will be shown below that 



j^P^ (iIpIU- + II (/) JL.) + II (/) Jl^. < a 



(5.3) 



9it)= (-.t) 



L'i(f7) 



as t ^ 00, 



(5.4) 



where 



/(t) 



1 



p^{x, t)dx. 



Now, we assume that (15.41) holds, and continue the proof of Proposition 15.11 In fact, 
the inequality (15. 3p and the Poincare-Sobolev inequality imply that 



c{n) 



< C 



p'-p'){:t) 



< c 



p'-p'){:t) 



2/3 

L'i{Q) 
2/3 



\\{p\{-M 



1/3 

L2 



0, as t — s> 00. 



(5.5) 



This suffices to finish the proof of Proposition 15.11 due to the following simple fact 

~f?{t) > -p\t) = p^^ = 1, for any t > 0. 

It remains to prove (15.40 . First, it follows directly from (15. 3p and the Poincare- 
Sobolev inequality that 



g(t)dt < C sup 

0<t<T 









Jo 



T 



\\{pW\hdt<c. 



(5.6) 



Next, we prove that 



\g'{t)\dt < C. 



(5.7) 
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Note that (12.91) as well as the boundary condition fl2.6p or (12 .yp imply that 



g'it) 



4(pM / {p'-p^Ydx 



-46 



y/p^udx - 4 ( j J (p- p'>) dx 



= h + h. 
It follows from Q and (l233ll that 

T pT 
\h\dt < C 







< C 



< c 



p'-7) {p\p'-''^^pudx 



dt 
dt 



p'-p' 



dt 



dt 



T 



< c / UpXWUt 

Jo 

< C. 



The uniform entropy estimate ( 12. 13^ . together with (12. 9p . gives that 



sup 

o<t<r 



h sup \pt)| 

0<t<T 



6 sup 

0<t<T 



[p''-^l^)J~''^pudx 

< c sup (||(p"''/'),ILJIv^«I|l^I 



< C sup 

0<t<T 



b—a 



< c. 



This together with Poincare inequality and (15.31) yields 



\I2\dt < C 



T 



P'-P' 



dt<C \M)x\\i2dt < C. 

Jo 



(5. 



(5.9) 



(5.10) 



(5.11) 



The estimate (15.71) thus follows directly from (I5.8I) - (I5.11I) . Hence the desired estimate 
(15. 4p follows from (15.61) and (15. 7p . The proof of Proposition 15.11 is completed. □ 
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5.2 Regularity and asymptotics of weak solutions for large 
time 

It is usually difficult to get information about the velocity field for the global entropy 
weak solution (p, ^/pu), in the sense of Definition [23] to the IB VP for the Compressible 
Navier-Stokes equations fl2.1l) - fl2.2p with initial data fl2.3l) and boundary values fl2.6p 
or (12.71) in the appearance of vacuum states. After vacuum states vanish, however, it 
will be shown that the velocity field u can be defined with enough regularity and the 
nonlinear diffusion term is represented in terms of the velocity u and the density p. 
The momentum equation becomes a uniform parabolic equation, and the weak solution 
(P) \/P'") ~ (P' \fP ' ^) becomes a strong solution. 

Proposition 15.11 implies that there is a time Tq > after which the density of the 
global entropy weak solution (p, u) to the IB VP problem for fl2.ip - fl2.3l) together with 
(12. 6p or (12. 7p is strictly positive and (p, u) satisfies the finite entropy estimate (I2.13p . 
Consider the IBVP problem (I2.ip - (l2.2p again for time t >Tq with data given at time 
t = To by 

p{x, To) = ^lim p{x, t), u{x, Tq) = ^hm ^==, (5.1) 

and note here that away from vacuum the Dirichlet boundary condition (12. 6p reduces 
to 

n(0,t) = n(l,t) = 0, t>To. (5.2) 

We then have the regularity property of the solution for the compressible Navier- 
Stokes equations (I2.ip - (l2.3p with the Dirichlet boundary condition (12. 6p or the periodic 
boundary condition (12.70 for large time. 



Proposition 5.2 Under the assumptions of Theorem \2.!% let (p, ^/pu) be the global 
entropy weak solution to the IBVP for the compressible Navier-Stokes equations i\2.1\i - 
i\2.2^ with initial data i\2.3\) and boundary value i\2.6^ or i \2.7 ) in the sense of Defini- 



tion \2.5l Let Tq > so that the global weak solution (p, u) satisfies for two positive 
constants p± that 

< p- < pix,t) < p+, V (x,t) G X [To,oo). (5.3) 



Then, (p, y^u) = (p, ^/p■u) is the unique strong solution to the IBVP for the compress- 
ible Navier-Stokes equations (\2. lf\ - (\2.2tl and (I5.jp with the boundary condition (\2.6^ 
or i \2. 7p for t > Tq. Moreover, the regularity (\2. 36^ and the long time behavior { \2.38^ 
hold. 



Proof: We only prove Proposition 15.21 for the Dirichlet case below, the periodic 
case can be treated similarly. 



^Note here that away from vacuum the Dirichlet boundary condition (|2.6p reduces to the usual 
one u(0, t) = t) = 0, t> Tq. 
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Step 1. Regularity. It follows easily from proposition 15.11 that there exist some Tq 
and a constant p_ > such that for all t >Tq, 

mlp{x,t) > p_ > (5.4) 

which, together with fl2.13p . implies 

peL'^{To,T;H\Q)) (5.5) 

for any T > Tq. By the continuity of p G (7(1] x [0, oo)) there exists some a > small 
enough such that for any t > Tq — a, 

inf p(x,t) > ^ > 0. 

This implies that one can define the velocity u for any global entropy weak solution in 
the sense of Definition 12.51 after the vanishing of vacuum states by 

u =: , for t >Tq — a. 

It then follows from the definition and (12.81) that 

M e L°°(ro -a,r;L2(l])) (5.6) 

for any T > Tq. Noting that fl232l) implies that for any ^{x) G C^ip^), tp{t) G Co"^(To- 
fT,T), 

ij{t) [ Ap-^fpdxdt 

T 



To~cr 

T 



^{t) I p"~'^^VP<p^''^)xdxdt 

a 

2a 



^{t) I ip''-'/^),^up-''ipdxdt 

To-a 



2a -1 

i/j{t) / ^up'^^^ifixdxdt 
To-cr Jn 

T f 

ip(t) / uip^dxdt, 
To- a Jn 

we can define the spatial derivative of velocity and, together with (12.130 . its regularity 

as 

u, = ^eL\nx{To-a,T)). (5.7) 

In terms of (15.61) . (15. 7p and (12.91) we are also able to justify the Dirichlet boundary 
condition (12. 6p for the velocity u 

u{0, t) = u{l, t) = 0, for any t > Tq - a. (5.8) 
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Thus, dES]), (Q and ([EHD show 

u G L^{To - a, T- H^{n)) n L°°(To - cr, T; ^^(fi)) (5.9) 

for the case of the Dirichlet boundary conditions. Note here that u G L'^{Tq — 
a, T; ifppj.(f2)) fl L°°(To — a, T; L'^{^)) in the case of periodic boundary conditions. We 
thus obtain from (12. 9p and (12.101) that the solution (p, u) satisfies 

Pt + (pM)x = a.e. in Qx (Tq -a,T), (5.10) 

and 



puiftdxdt + I I i^pu"^ — p^'Ux + p'^) ^xdxdt = (5-11) 

for any ip{x,t) G C^{Q x (Tq — cr, T)) for the Dirichlet case. The Eq. (15.111) can be 
re- written in terms of (15.91) . (I5.10p and (I2.13P as follows 

/ / {ucpt - p"~^Ux^x + (p"~^Px - u) Uxip)dxdt = / 'yp^~'^pxipdxdt (5.12) 

Jto-ct Jn JTa-cr Jn 

for any (p{x,t) G C^{fl x (Tq — o", T)) for the Dirichlet case. 

Noticing that p^'^p^ -uE L°°(To - a,T; L'^{n)) due to ([53D and and 

using standard regularity results for linear parabolic equations (see [25]), we get that 

u G L2(To, T; i/2(fi))) n H\To, T; L^Vl)) (5.13) 

for the Dirichlet case. It is noted here that u G L'^{To,T] H^^^{Q)) n H\To,T] L'^ift)) 
for the periodic case. 

Step 2. Uniqueness. We shall show that if there exists another solution {rj, v) to the 
compressible Navier-Stokes equations (12. II) -( 12^ with the following initial data and 
Dirichlet boundary conditions 



{r],v){x,To) = (p,u)(x,To), 
v{0,t) = v{l,t) = 0, 



(5.14) 



such that 



p^<r]eL^{To,T;H\Q)), 

V G L-(To, T; H\n)) n L\To, T; H^Q)) n H\To, T; L^Q)), 



(5.15) 



then 

p = rj, u = v a.e. in Q x {Tq,T). (5.16) 
In fact, it follows from fl5A0ll . flHTTl) and (ICTl) that 

- j pu'^dx + I I p'^uldxds — I I p^Uxdxds = - f pu^{x,TQ)dx (5.17) 
2 Jn Jto Jn Jto Jn 2 Jq 
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for all te{To,T), while fISTTD and flSTTSD imply that 

puvdx + / / p°'UxVxdxds — / p^'v^dxds 
n Jto Jn Jto Jn 

puv{x,TQ)dx+ / / pu{vt + uVx)dxds (5.18) 
n Jto Jn 

for all t G (To,T). To estimate the second term on the right hand side of flS.lSp . we 
use the decomposition 

pvt + puv^ = {p- V){'"t + vv^) + p{u - v)v^ + {'q°'v^)^ - {t]'^)^. (5.19) 

Multiplying 05.191) by u and we integrating by parts give 

/ pu{vt + uVx)dxds 
To Jn 

{p — f])u{vt + vvx)dxds + / pu{u — v)vxdxds 
To Jn Jto Jn 

— / rj°'vxUxdxds + / rfuxdxds. (5.20) 
Jto Jn Jto Jn 

Substituting into fl5JB gives, for a.e. t E (To,T), that 

puvdx + / p°'UxVxdxds — / {p^Vx + rfux)dxds 
n Jto Jn J To Jn 

puv{x, Tojdx + / (p — f])'^{'"t + vVx)dxds 
n Jto Jn 

/ pu{u — v)vxdxds — / ■r]°'VxUxdxds. (5-21) 
To Jn Jto Jn 

Multiplying (15.191) by v and integrating the result over fl x (TQ,t) lead to 

- j pv'^dx = - I pv^{x,TQ)dx + I I {p — ri)v{vt + vVx)dxds 
2 Jn J To Jn 

+ pv{u — v)vxdxds — I I rj"'vldxds 

Jto Jn J To Jn 



rpVxdxds. (5.22) 



iTo Jn 

We obtain after adding up fl5.17p and 05.221) and subtracting 05.2ip that for all 
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te{To,T), 



1 

2 



p{u — vYdx+ I / p°'{u — v)l.dxds 
I J To Jn 

= / {p" - r]°'){v - u)xVxdxds + / / {p - 7]){vt + vvx){v - u)dxds 
Jto Jn Jto Jn 

~ / ~ vYv^dxds —I I {p^ — rf) {v — u)xdxds 
J To Jn J To Jn 

< Ce f (Hp" - tWUKWl- + Hp - r^WlAW + w^wl, + \\p^- r^^\\l,) ds 

J To 

+C I ||fx||L°° / p{u — vY dxds + C e I \\u — v\\\ids. (5.23) 
Jto Jn Jto 

Next, we estimate the term 

9{t) = \\p-v\\h + W-t\\l^ + \\p'-V'\\l^- 
For any /3 > 0, fl5J[5|) . (I57[3|) and imply that 

- + vifJ" - + {u- v){p% + PpP{u - v), + (3{p^ - ri^)v, = 0. (5.24) 
One can derive from this that 

iWp^-v'Wh), < IIP^-/lli.(c|K^.||L^ + c,||p.||i. + c,) 

+e\\u — v\\1oo + e\\{u — f)z||L2. (5.25) 
Thus, since 5'(To) = 0, we obtain from (15.251) with /5 = l,a,7 respectively that 

9{t)< [ 9{s){C\\v4L^+CMl, + C,)ds + Ce [ \\u-v\\l,ds. (5.26) 
Jto JTq 

Now fl516D is a consequence from fl5:23|) . ([SJSD, fl535D . fl53|) and ([EID- The proof 
of large time convergence (I2.38P follows directly from the standard arguments (see [44] 
for instance) with the help of entropy inequality (I2.13p . The proof of Proposition 15.21 
is completed. □ 



5.3 Finite time blow-up 

In this subsection, we shall prove the Theorem 12.41 about the finite time blow- 
up phenomena as an immediate consequence of Theorem 12.31 Propositions 15.11 and 
Proposition 15. 2[ 
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Proof of Theorem\2.4\ We will prove (12.421) only. The proof of (I2.43P is similar. If 
(I2.42P fails, then there exists a fixed constant ?7 > 0, such that 

rTi+ri 

/ ||iix||L°°'^s < oo. (5-27) 

For any {x,t) E ft x (Ti,Ti + 77], the particle path x{s) = X(s;t,x) through (x, t) is 
given by 

'^X{s;t,x)=u{X{s;t,x),s), < s < t < T, + 7], ^ ^g) 
X{t; t,x) = X, Ti < t < Ti + 7], X E il, 

which is well-defined due to (I5.27P and (I2.36p . Consequently, one obtains via a standard 
argument from the transport equation (12. ip 

p(a;,t) = p(X(Ti;t,x),Ti)exp |- y Uy{y, s)\y=xis;t,x)ds^ (5.29) 

for any {x,t) E Q x (Ti,Ti + rj]. On the other hand, it follows from (I5.27P and (I5.28P 
that for any x E Q there exists a trajectory x = x{t) E Q for t E [Ti,Ti + 77] so 
that X{Ti]t,x{t)) = X. In particular, there exists a trajectory x = Xi{t) E Q for 
t E [Ti,Ti + T]] so that X{Ti;t,x(t)) = Xi with (xi,Ti) determined by (I2.40p . namely, 
p(xi,Ti) = 0. Thus, due to fl5.29p . we deduce from f l5.27p that 

p{xi{t),t) = for all t E (Ti, Ti + r/], 

which contradicts (12.401) . Thus, the blowup phenomena (I2.42p happens. The proof of 
the Theorem 12.41 is completed. □ 
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